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Abstract. 

In the framework of the theory of open systems based on completely positive 
quantum dynamical semigroups, we give a description of continuous variable Gaussian 
Renyi-2 (GR2) quantum discord for a system consisting of two non-interacting non- 
resonant bosonic modes embedded in a thermal environment. In the case of both an 
entangled initial squeezed vacuum state and squeezed thermal state, GR2 quantum 
discord has non- negative values for all temperatures of the environment. We describe 
also the time evolution of GR2 classical correlations and GR2 quantum mutual 
information of the quantum system. All these quantities tend asymptotically to zero 
in the long time regime under the effect of the thermal bath. 



1. Introduction 

In recent years there is an increasing interest in using non-classical entangled states 
of continuous variable systems in applications of quantum information processing, 
communication and computation [UI21E]- In this respect, Gaussian states, in particular 
two-mode Gaussian states, play a key role since they can be easily created and controlled 
experimentally. Due to the unavoidable interaction with the environment, in order to 
describe realistically quantum information processes it is necessary to take decoherence 
and dissipation into consideration. Decoherence and dynamics of quantum entanglement 
in continuous variable open systems have been intensively studied in the last years 



In quantum information theory an interesting family of additive entropies is 
represented by Renyi-a entropies [21 J, defined by 



In the limit a — > 1 they reduce to the von Neumann entropy S(p) = — Tr(plnp), 
which in quantum information theory quantifies the degree of information contained 
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in a quantum state p, and in the case a = 2 we obtain ^(p) = — lnTr(p 2 ), which is 
a quantity directly related to the purity of the state. Von Neumann entropy satisfies 
the strong subadditivity inequality which is a key requirement for quantum information 
theory, implying in particular that the mutual information, which measures the total 
correlations between two subsystems in a bipartite state, is always non-negative. Renyi 
entropies are powerful quantities for studying quantum correlations in multipartite states 
[2"2"[ |2"B"] . In general Renyi-a entropies for a ^ 1 are not subadditive. In Ref. [21] it 
was demonstrated that Renyi-2 entropy provides a natural measure of information for 
any multimode Gaussian state of quantum harmonic systems. It was proven that for all 
Gaussian states this entropy satisfies the strong subadditivity inequality, which made 
possible to define measures of Gaussian Renyi-2 (GR2) entanglement, total, classical, 
and discord-like quantum correlations based on this entropy. In this sense one could 
regard Renyi-2 entropy as a specially meaningful choice to develop a Gaussian theory 
of quantum information and correlations |24j. 

In this paper we study, in the framework of the theory of open systems based on 
completely positive quantum dynamical semigroups, the dynamics of continuous variable 
GR2 quantum discord of a subsystem consisting of two bosonic modes (harmonic 
oscillators) interacting with a common thermal environment. In Sec. 2 we give the 
general solution of the evolution equation for the covariance matrix, i.e. we derive the 
variances and covariances of coordinates and momenta. We are interested in discussing 
the correlation effect of the environment, therefore we assume that the two modes are 
independent, i.e. they do not interact directly. The initial state of the subsystem is taken 
of Gaussian form and the evolution under the quantum dynamical semigroup assures the 
preservation in time of the Gaussian form of the state. In Sec. 3 we show that in the case 
of both an entangled initial squeezed vacuum state and squeezed thermal state, GR2 
quantum discord has non-negative values in time for all values of the temperature of the 
environment. We describe also the time evolution of GR2 classical correlations and GR2 
quantum mutual information, which measures the total amount of GR2 correlations of 
the quantum system.. All these quantities tend asymptotically to zero in time under 
the effect of the thermal bath. A summary is given in Sec. 4. 

2. Equations of motion for two modes interacting with an environment 

We study the dynamics of a subsystem composed of two non-interacting bosonic modes 
in weak interaction with a thermal environment. In the axiomatic formalism based on 
completely positive quantum dynamical semigroups, the Markovian irreversible time 
evolution of an open system is described by the Kossakowski-Lindblad master equation 
[251 [26J. We are interested in the set of Gaussian states, therefore we introduce such 
quantum dynamical semigroups that preserve this set during time evolution of the 
system. The Hamiltonian of the two uncoupled non-resonant harmonic oscillators of 
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identical mass m and frequencies oj% and U2 is 
H 



1 (pI+pD + j^ + uIv 2 ), 



(2) 



2m xrx ' ry 

where x, y are the coordinates and p x ,Py are the momenta of the oscillators. 

The equations of motion for the quantum correlations of the canonical observables 
x, y and p x ,p y are the following (T denotes the transposed matrix) 
da(t) 



dt 



Ya{t) + a{t)Y l + 2D, 



(3) 
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and the dissipation constant A are real 
quantities. We introduced the following 4x4 bimodal covariance matrix: 



a(t) 



( o- xx {t) o- xpx {t) o- xy (t) 

(*) °jtafl„(*) Via*® 
(t) 



a 



X P y 



a 



XP.T 



a 



xy 



(t) 



Px 

T yp 

'p x p y 



a 



<T„ Vx (t) 



a 



yy 
ypy 



PxPy 

ypy 



(t)\ 
(t) 



<T.. m „,(t) 



(Tun., (t) 



a 




(5) 



PyPy 



(t)J 



where A, B and C are 2x2 Hermitian matrices. A and B denote the symmetric 
covariance matrices for the individual reduced one-mode states, while the matrix C 
contains the cross-correlations between modes. The elements of the covariance matrix 
are defined as = Tr[p{RiRj + RjRi}]/2,i, j = 1, ..,4, with R = {x,p x ,y,p y }, which 
up to local displacements fully characterize any Gaussian state of a bipartite system [27J . 
All the measures defined below are invariant under local unitaries, so we will assume 
our states to have zero first moments, (R) = 0, without loss of generality. 
The time-dependent solution of Eq. ^ is given by [26] 

a(t) = M(t)[a(0) - a(oo)]M T (t) + a(oo), (6) 

where the matrix M(t) = exp(Yt) has to fulfill the condition lim^oo M(t) = 0. The 
values at infinity are obtained from the equation 

Ya(oo) + a(oo)F T = -2D. (7) 



3. Dynamics of GR2 correlations 

3.1. Time evolution of GR2 quantum discord 

The GR2 mutual information T 2 for an arbitrary bipartite Gaussian state pab is defined 
by 

^h(Paib) = S 2 {p A ) + S 2 {p B ) - <S 2 (pab), (8) 
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where pa and ps are the two marginals of pab- It was shown [21] that a-.b) > and 
it measures the total quadrature correlations of pab- 

One can also define a GR2 measure of one-way classical correlations Ji{pa\b) 
[2H [28] , as the maximum decrease in the Renyi-2 entropy of subsystem A, given a 
Gaussian measurement has been performed on subsystem B, where the maximization 
is over all Gaussian measurements, which map Gaussian states into Gaussian states 



Following [31] and the analysis of Gaussian quantum discord using von Neumann 
entropy [29, 30J, in [24J it was defined a Gaussian non-negative measure of quantumness 
of correlations based on Renyi-2 entropy, namely GR2 discord, as the difference between 
mutual information (|8l) and classical correlations: 



Up 



A:B 



B 



(9) 



^2 (pA\B 

For a general two- mode Gaussian state pab, with A and B single modes, closed formulae 
have been obtained for GR2 classical correlations and GR2 discord [24J. For pure 
bipartite Gaussian states pab the following equalities hold: ^^(pa-.b) = 3i{.Pa\b) = 
^>2{pa\b) = S 2 (pa)- 

The covariance matrix a ^ of any two-mode Gaussian state can be transformed, 
by means of local unitary operations, into a standard form of the type [32] 
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where a, b > 1 



(10) 



> 0, and we can set 



(a 2 - 1) (b 2 -l)-ab (c 2 + c\) - 2c_c+ + c 2 _c\ 
c + > |c_| without losing generality. These conditions ensure that the uncertainty 
relation a > zw e2 /16 is satisfied (cj = is the symplectic matrix), which is a 

requirement for the covariance matrix a to be associated with a physical Gaussian state 
in a two-mode infinite-dimensional Hilbert space [32J. For pure Gaussian states, b = a, 
c + = — c_ = \J a 2 — 1. 

For generally mixed two-mode Gaussian states Pab, the Renyi-2 measure of one- 
way quantum discord ^ has the following expression if the covariance matrix a is in 
standard form (Hoi El: 



T>2{pa\b) 



In b- -ln(detC) + ^lne 2 , 

_ 



(11) 



where 
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a I a — f 



, if (ab 



2 c 2 



a + bc 2 _)) (ab 2 c\ - c 2 (a + bc 2 + )) < 0; 



2\c-c + \J(a{b 2 -l)-bc 2 _)(a(b 2 -l)-bc 2 + ) + (a(b 2 ~l)~-bci )(a(b 2 -l)-bc 2 + )+c 2 _c 



These formulae, written explicitly for standard form covariance matrices, can be 
reobtained in a locally invariant form by expressing them in terms of the four local 



otherwise. 



symplectic invariants of a generic two-mode Gaussian state 



h = 4detA, I 2 
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4det.B, ^3 = 4detC, J4 = 16 deter [34] . This can be realized by inverting the relations 
Ii = a 2 ,I 2 = b 2 ,I 3 = c + c_,Ii = (ab — c + )(ab — c_) [32]. The obtained expressions 
are then valid for two-mode covariance matrices in any symplectic basis, beyond the 
standard form. 

We assume that the initial Gaussian state is a two-mode squeezed thermal state, 
with the covariance matrix of the form [351 ES] 

/ a s c s \ 



*.(0) = \ 



a s — c s 
c s b s 



(13) 



V -c a b s J 
with the matrix elements given by 

a s = n\ cosh 2 r + n 2 sinh 2 r + - cosh 2r, (14) 
b s = ni sinh 2 r + n 2 cosh 2 r + - cosh 2r, (15) 

c s = + n 2 + 1) sinh2r, (16) 

where ni, n 2 are the average number of thermal photons associated with the two modes 
and r denotes the squeezing parameter. In the particular case m — and n 2 = 0, 
(13) becomes the covariance matrix of the two-mode squeezed vacuum state [37J. A 
two-mode squeezed thermal state is entangled when the squeezing parameter r satisfies 
the inequality r > r s [35], where 

cosh 2 r,= ( " 1 + 1)( " 2 + 1) . (17) 
^1 + ^2 + 1 

We suppose that the asymptotic state of the considered open system is a Gibbs state 
corresponding to two independent bosonic modes in thermal equilibrium at temperature 
T [2H]. Then the quantum diffusion coefficients have the following form (we put h = 1) 

muxD xx = = - coth — — , 

moj\ 2 2kl 

'^ D » = ^ = r oth ^ < 18) 

D X p x — Dyp y = D X y = Dp xPy = D'XPy = D yp^ = 0. 

The evolution of the GR2 quantum discord is illustrated in Figs. 1 and 2, where 
we represent the dependence of T>2(pa\b) on time t and temperature T for an initial 
entangled Gaussian state with the covariance matrix given by Eq. (13), taken of the 
form of a two-mode squeezed vacuum state, respectively squeezed thermal state, and for 
such values of the parameters that satisfy for all times the first condition in formula (12). 
The GR2 discord has positive values for all finite times and in the limit of infinite time 
it tends asymptotically to zero, corresponding to the thermal product (separable) state, 
with no correlations at all. We also notice that the decay of GR2 discord is stronger 
when the temperature T and dissipation constant A are increasing. 
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Figure 1. GR2 quantum discord T>2 V {pa\b) versus time t and temperature T for 
an entangled initial squeezed vacuum state with squeezing parameter r — 0.23, 
rii = 0,n 2 = and A = 0.1, u>i = 1,cj 2 = 2. We take m = % = k = 1. 




Figure 2. GR2 quantum discord T>2(pa\b) versus time t and temperature T for 
an entangled initial non-symmetric squeezed thermal state with squeezing parameter 
r = 0.1, n\ = 1, 712 = and A = 0.1, wi = 1, u)2 = 2. We take m = h = k = 1. 



5.1?. GR2 classical correlations and quantum mutual information 

For a mixed two-mode Gaussian state Pab, GR2 measure of one-way classical 
correlations J2{pa\b) has the following expression if the covariance matrix a is in 
standard form ([l0[ gl| : 

Ji{.Pa\b) = hia - - hi£ 2 , (19) 
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Figure 3. GR2 classical correlations Ji{_Pa\b) versus time t and temperature T for 
an entangled initial non-symmetric squeezed thermal state with squeezing parameter 
r = 0.1, ni = l,7i2 = and A = 0.1, wi = l,ui2 = 2. We take m = h = k = 1. 




Figure 4. GR2 quantum mutual information T-i^pa-.b) versus time t and temperature 
T for an entangled initial non-symmetric squeezed thermal state with squeezing 
parameter r — 0.1, n\ = 1, n2 = and A = 0.1, = 1, CJ2 = 2. We take m = K = k = 1. 



where £2 is given by Eq. (12), while the expression of the GR2 quantum mutual 
information, which measures the total correlations, is given by 



1i{pa:b) = In a + In b - - m(det C). 



(20) 



In Figs. 3 and 4 we illustrate the evolution of classical correlations J%{pa\b) 
and, respectively, quantum mutual information X 2 (pa : b) as functions of time t and 
temperature T for an entangled initial Gaussian state, taken of the form of a two- 
mode squeezed thermal state (13), and for such values of the parameters that satisfy 
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Figure 5. GR2 quantum mutual information, classical correlations and quantum 
discord versus time t and temperature T for an entangled initial non-symmetric 
squeezed thermal state with squeezing parameter r = 3, ri\ = 3, tt-2 = 1 and 
A = 0.1, Mi = l,u)2 = 2. We take m = K = k = 1. 



for all times the first condition in formula (12). Both these quantities manifest a 
similar qualitative behaviour similar: they have positive values for all finite times and 
in the limit of infinite time they tend asymptotically to zero, corresponding to the 
thermal product (separable) state, with no correlations at all. One can also see that 
the classical correlations and quantum mutual information decrease with increasing the 
temperature of the thermal bath and the dissipation coefficient. For comparison, GR2 
mutual information, classical correlations and discord are represented also on the same 
graphic in Fig. 5. In the considered case the value of GR2 classical correlations is larger 
than that of quantum correlations, represented by the GR2 quantum discord. 



4. Summary 

We investigated the Markovian dynamics of Gaussian Renyi-2 quantum correlations for a 
subsystem composed of two non-interacting bosonic modes embedded in a thermal bath. 
We have analyzed the influence of the environment on the dynamics of GR2 quantum 
discord, classical correlations and quantum mutual information, which measures the 
total GR2 correlations of the quantum system, for initial squeezed vacuum states 
and non-symmetric squeezed thermal states, for the case when the asymptotic state 
of the considered open system is a Gibbs state corresponding to two independent 
quantum harmonic oscillators in thermal equilibrium. The dynamics of these quantities 
strongly depend on the initial states and the parameters characterizing the environment 
(dissipation coefficient and temperature). Their values decrease asymptotically in time 
with increasing the temperature and dissipation. 

The study of time evolution of these measures may contribute to understanding the 
quantification of correlations defined on the base of the Renyi-2 entropy, which could 
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represent a special tool in developing a Gaussian theory of quantum information and 
correlations. 
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